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ANALYSIS  OF  OEKSRAL  DYNAMIC  3YSTEKS 
IHCUS?X!3  OHS  Nr^LIESAR  5LEKSKT 

Introduction  _^A_D^5Si^ 

la  a previous  report*1-  describing  the  appll  Ci>wlon  o*  block 
diagram  methods2"^  to  nonlinear  circuit  problems  , it  was  shown 
that  the  mathematical  relations  for  a variety  of  circuits  in- 
cluding one  nonlinear  element  could  be  indicated  by  a block 
diagram  of  the  type  given  in  Fig.  1.,  In  this  s-  O^Cs), 

Q2(s}f  and  fi  i ; , re  transfer  functions^  which  express  in 
operational  form  dynamic  relations  which  could  also  b®  given 
as  linear  differential  equations*  The  forward  block  represents 


INPUT  BLOCK  FORWARD  BLOCK  OUTPUT  BLOCK 


Figure  1.  General  Block  Diagram 


a function,  f(x),  which  describes  an  instantaneous  nonlinear 
relation  between  the  variables  x and  y,  It  was  also  pointed 
out  that  circuits  exist  whose  block  diagrams  cannot  be  reduces 
to  this  general  form. 

Further  study,  euxamarized  in  this  report,  indicates  that 
networks  of  a fairly  general  class  can  be  described  in  this  way. 


i 

% 


_ 2 


A theorem  13  presented  which  permits  the  I'c-rw  arc  f<nd  feedback 
tlvCKS  to  bo  Interchanged,  if  this  should  prove  desirable*  The 
exceptions  are  likewise  found  to  constitute  a general  class  to 
which  similar  methods,  with  some  complications.,  - ~n  be  applied. 

As  Indicated  previously-*,  the  variable  bo  computed 

from  e by  conventional  methods?  If  y can  be  f-.-rnd,  eQ  can  also 
be  computed  as  the  response  of  a linear  system.  The  real  problem, 
therefore,  is  the  analysis  of  feedback  system  ith  a nonlinear 
forward  block  and  a linear  feedback  block,  or  \ v*irsa*  In 
this?  report,  a possible  str-p-by-atep  method  for  t his  purpose  la 
described.  *fhe  method  involves  the  dsnsrminatiso*^  of  by  oper- 
ational methods,  solution  of  the  feedback  system  by  a combination 
of  an  approximation  to  the  superposition  integrr".  and  a graphical 
solution  of  two  simultaneous  algebraic  equation*.  suod  finally  & 
calculation  of  eQ  by  any  appropriate  method.  A f~r  examples  are 
included  to  illustrate  the  method  and  the  quest 5 u*.  of  accuracy 
is  considered  briefly. 

Although  the  examples  discussed  are  problems  in  electric 


circuit* 


control  systems , the  procedures  are  applicable  to 


other  types  of  dynamic  systems  by  the  well-known  Loounique  of 
electro -mechanical  analogies  or  by  direct  construction  of  the 
block  diagrams. 


VnsmNB 


The  block  diagrams  for  individual  linear  elements  are  given 
in  Fig.  2.  Two  sets  of  diagram  a , called  "impedance"  and  *’ admitt- 
ance" , are  required  because  of  the  possibility  that  either 


voltages  or  currents  may  appear  as  Inputs  to  the  blocks , Mutual 
inductance  and  the  amplification  factor  or  transconduct  ones  of 
vacuum  tubee  may  also  be  required  in  the  block  diagrams* 


SSjEHEN 
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Inductance 


Capacitance 


Figure  2.  Block  Representation  of  Linear  Elements 


She  block  diagram  representation  of  nonlinear  elements  la 
not  so  simple j since  it  generally  requires  the  specification  of 
a nonlinear  function  and  a dynamic  relation.  'Hie  representations 
adopted  for  the  common  circuit  elements  are  given  in  FAg6  3*  2be 
script  symbols  \jR3%,&9  end  sc  forth)  denote  functions,  with 
the  output  of  the  block  ^nsidered  as  a faction  of  the 

Input.  This  function  ehjt  he  given  analytically  cr  may  only  be 
expressible  in  graphic ai  or  tabular  form. 

For  a nonlinear  resistor,  the  voltage -current  relation  may 
be  given  as 

s =(&(i)  (1* 

or  i ==  (2) 

depending  on  whether  the  current  or  voltage  is  considered  to  bo 
the  independent  variable  or  input. 
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A nor.lir.aar  inductor  may  be  characterised  1 sy  two  relations: 


1 

P 


§il). 


(3) 


which  represents  the  saturation  or  magnetisation  curve , and 

& 

9 ” ^ dt  * 

which  ia  Faraday * s lajr*  If  several  windings  appear  on  the  sami? 
cor#,  the  flux  is  taken  as  a function  of  the  net  ampere  -turns  * 
A nonlinear  capacitor  is  also  described,  top  two  relations: 


(4) 


Q<®>, 


q = WO).  (5i 

which  is  similar  to  the  saturation  curve  of  a magnetic  material’, 
and 

(6) 


V 


I 


a basic  definition* 

It  may  be  noted  in  passing  that  the  principle  of  duality  is 
applicable  in  nonlinear  circuits  and  that  the  saturation  curves 
©f  ferromagnetic  and  ferroelectric  materials  correspond  closely 
enough  that  solutions  for  one  problem  can  sometimes  be  used 
directly  as  solutions  of  the  associated  dual  network. 

PiagpsB 

The  block  diagram  for  any  network  having  the  configuration 
shown  in  Fig.  4(a)  or  5(a)  can  be  reduced  to  the  general  form 
given  in  Fig.  1.  In  those  networks  f the  nonlinear  element  appears 
between  two  four-terminal  networks,  in  either  a series  or  shunt 
position.  Although  nonlinear  resistors  are  shown,  any  of  the 
other  nonlinear  elements  may  be  substituted. 
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Hi©  four-tersinaX  networks  may  be  characterized  by  their 

equivalent  T-  or*  ?r-aections , whichever  is  most  convenient-  as 

shown  in  Figs.  4(b)  and  5(b).  The  impedances  of  tbs  generator 

supplying  the  network  and  the  measuring  device  at  the  outFdt 

are  assumed  to  be  included  in  the  networks  N and  H*  Should  it 

be 

be  desirable*  the  networks  can  aXso/speeifted  by  their  open- 
circuit  Impedance  parameter  cr  the  short-circuit  admittance 

•7 

parameter's  . 

The  original  block  diagrams  , based  on  the  circuit  equations  * 
are  given  In  Figa.  4(c)  snd  5(c)*  Th©  final  diagrams,  obtained 

•>  3. 

by  a few  simple  block  transfonaations*1"^,  are  given  in  Figs.  4(d) 
and  5(d).  Study  of  those  final  diagrams  indicates  that  they  can 
also  be  obtained  by  direct  application  ef  standard  network  theorems; 
Fig.  4(d)  can  be  derived  from  the  circuit  diagram  by  the  use  of 
Thavsnin's  theorem,  and  Fig.  5(d)  by  Horton's  theorem. 

It  will  be  observed  that  a series  nonlinear  element  is  beat 
described  by  its  admittance  function  (voltage  input  and  current 
output),  and  that  a shunt  element  is  described  by  an  impedance 
relation.  Likewise,  it  nay  be  seen  that  4^  does  not  appear  In 
Fig*  4(d)  and  YA  does  not  appear  in  Fig.  5(d). 

The  same  type  of  final  block  diagram  is  obtained  if  the 
four -terminal  networks  are  coupled  by  a nonlinear*  (iron-core) 
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If  the  series  element  is  shunted  or  bridged  by  an  unlike 

linear  element,  or  if  an  unlike  linear  element  is  placed  in 

series  with  the  shunt  element,  a network  is  produced  whose 

not 

block  diagram  does/ reduce  to  the  general  fora  cf  Figs  1.  oucb 


networks  are  shown  ir»  figs.  6 and  T. 


Figure  6.  Par&llal  Combination  between  Two  Networks 
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A possible  method  for  dealing  with  exceptions  of  thia  type 
is  suggested  In  Fig.  8,  A circuit  diagram  for  the  network  of  / 
Fig,  6f  using  T -section  represent  at  ions  for  the  four-terminal 
networks  and  expressing  the  elements  in  terns  of  the  impedance 
parameters 3 is  gi yen  in  Fig*  61a)*  In  developing  the*  block 
diagram  given  in  Fig,  3(b)  * Thevenir*s  theorem  is  applied  to 
the  left-hand  network  ( N 3 to  obtain  an  open-circuit  voltage, 
e(Zj2/2n> , and  an  equivalent  source  .impedance , l/Xgo  * l;ke 
voltage  across  the  parallel  combination  of  Yp  end  the  nonlinear 
element  is  th6n  computed  by  subtracting  the  voltages  across 
( 1/^22  i ^nd  from  the  open-circuit  voltage.  rfith  the  voltage 

across  the  parallel  combination  available*  the  currant  entering 
tbo  rrght -hend  network  (N)  can  be  computed  and  then  used  to  find 
the  outjwt.  voltage,  as  well  as  the  internal  voltage  drops. 

For  convenience,  the  feedback  block  in  Fig,  8(b)  is  c ailed 
w p"  in  developing  the  final  diagram  given  in  Fig,  8(c).  This 
diagram  can  be  obtained  by  block  transformations  or  by  applying 
Horton's  theorem  to  the  equivalent  circuit  from  which  Fig.-,  8(b) 
was  derived.  The  final  diagram  bears  some  resemblance  t-e  Fig,  1, 
although  the  nonlinear  block  row  appears  in  the  feedback  path, 
and  can  be  bandied  by  similar  methods* 

The  network  shewn  in  Fig.  7 may  be  treated  ir.  essentially 
the  seme  way  and  the  details  are  not  included  here. 


The  block  t rsnsf  crmat ien  shewn  in  Fig.  $?(»),  although  not 
widely  known  or  used*  can  easily  he  checked  by  solving  the  basi© 
equations  of  the  feedback  systems  or  by  application  of  tlxe  usual 
relation 
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It  the  forward  block  is  nonlinear,  the  transformation  shown 
in  Fig,  9(b)  is  applicable 3 this  follows  from  the  fact  that  the 
equations  sisy  be  written  as 


or  as 


X 

“ 

X1  ' *2 

(8) 

y 

= 

fix) 

(9) 

*2 

— i g 

jS(s)  y 

(10) 

*2 

= 

K1  * x 

(11) 

y 

^ ( s } 

(12) 

X 

— 

s(y) 

(13) 

where  the  function  g is  related  to  f in  such  a way  that 

x = g O)]  - x (14) 

It  m«y  be  said  that  the  forward  and  feedback  blocks  can  be 
interchanged  by  the  substitution  of  their  inverses  or  reciprocals* 
Application  of  this  theorem  permits  the  feedback  loop  in  Fig.  8(e) 
to  be  reduced  to  the  standard  fernu 

An  Illustration  of  the  theorem  is  presented  in  ^igs  10, 


where  a apes if ie  circuit  problem  is  considered. 


S. 


* 
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©w»  analytical  problem  pcsM  or  £%•  I hsg  ha  ststr-d  as 
fcJ.lowsj  given  a complete  defrey'^blen.  or  the  lisaa?  and  son.- 
linear  portions  of  the  ays  ter-  *"A  s.  ep^oifSsa  -inp^t,  *-;tf  * tf&ss- 
is  the  corroepcr^ir**  outp-  «.»  4 (t)?  A gssoral  app$o&*h  this 
probl*?**  is  described  sad  tilust  ***&*%&  in  ti«.  isasind&r'ef  the 
report. 

The  first  step  in  the  eolation  is  the  determination  of 
x,  (t) . which  may  be  accomplished  by  any  of  several  assthodp. 

Since  the  input  block  la  lir^  r,  the  differential  equation 
relating  e(t)  and  x^(t)  car  be  recovered  and  solved  by  classical 
methods . Using  operational  uatheds,  we  would  have 


xx(t)  = £~x  [&iCs)  E(s)j  (15) 

where  B(s)  is  the  transform  of  e(t).  If  the  steady-state  solution 
for  a periodic  input  la  required,  the  transfer  function  G^Cs) 
can  be  converted  into  an  expression  for  magnitude  and  phase 
at  reel  frequencies  by  the  substitution  of  Jo  for  s.  In  the 
general  case,  wh oro  neither  L&plc-ce  nor  Fourier  transforms  are 
especially  convenient,  it  is  possible  to  write® 


xx(t)  = f a£t  - X ) Bi(V  (16) 

*0 

where  g,(t)  = T G-„  (a)”j  „ the  impulse  response  or  weighting 

* «-  x - • 

function  for  the  input  block,  and  to  dot ermine  x1(t ) approxi- 
mately to  any  desired  accuracy. 


~ 16  - 


The  real  problem  is  the  calculation  of  the  response  of 


tils  system  ir  that  i«,  the  determination  or  y(t)  fee 


& given  x.s(t).  Using  Sqs* 
in  the  t.AS*-e  domain  a©  «*: 


(S)-(10)#  this  probles  may  be  written 
ir-tegra!  equation 


*(t) 

= X, (t)  - 

tt 

1 A (V  f 

[x(t  -<t)f  do 

(17) 

or 

x(t) 

= Xx(t)  - 

y 

f fid  -t) 

t ^x(tf|  &X 

(13) 

’-'O 

where 

/3(t) 

» a^Oujj, 

the  impulse 

response  or  weighting 

function  of  the  feedback  block*  Equations  (17)  end  (18)  sire  both 


Bpeeial  cases  of  Laleseo's  nonlinear  integral  equation^,  which 
has  been  applied  to  circuit  problems  by  Keller1**  and  Pipes11. 

In  the  analytic  treatment  of  this  type  of  equation,  a sequence 
of  functions  is  obtained  whose  limit  is  the  solution-  The  main 
difficulty  with  this  method  is  the  complexity  of  the  integrands 
encountered.  The  necessity  for  an  analytic  description  of  the 
nonlinearity  may  also  be  troublesome  * 

Approximate  solutions  of  the  system  equations  can  be  obtained 
by  step-by-step  methods,  using  either  one  of  two  procedures  based 
on  Eqs.  (8}-(10),  v/hich  will  probably  be  simpler  to  apply  than 
the  iteration  method  Just  described.  In  the  first  method,  the 
feedback  block  is  handled  by  writing  3q«  (10)  as 

r* 

x2(t)  = | /5(t)  y(t  - X ) d'fc  ; (19) 

^ b 

this  integral  is  then  evaluated  numerically  by  any  of  the  well- 
known  methods  for-  approximate  integration'**""- This  approach 
has  been  used  by  Mazelcky  and  Diederieh1^  in  wo- rk  on  strictly 
linear  systems. 


In  the  second  method;  y(t } is  resolved  into  component 
•signals  {step  functions,  r&ctsngul&r  or  triangular  pulses,  end 
e©  .forth)  for  which  the  response  of  the  feedback  block  can  be 
computed;  the  variable  x^Ct ) is  obtained  by  superposition  of 
the  responses  to  the  separate  components  of  y(t ) . This  method 
has  been  applied  to  nonlinear  systems  by  Tustin*' "J‘w,  Madwed-'9, 
end  H&rraer  v;  some  of  the  mathematical  questions  involved  have 
been  diacussod  by  Bro In  these  applications,  the  sequence 
of  numbers  representing  the  response  of  the  linear  portions  of 
the  system  to  the  basic  Input  component  has  been  found  by 
substituting  difference  operators  for  different ial  operators, 
rather  than  by  direct  calculation. 

The  first  method  will  be  employed  in  this  report,  since 
it  appears  to  offer  some  advantages.  Unlike  the  ether  methods, 
no  attempt  is  made  to  approximate  the  variables  in  components 
or  to  express  the  operation  of  differentiation  in  terms  of 
differences.  Errors  peculiar  to  the  first  method  are  made  in 
the  approximate  integration,  with  the  result  that  their  else 
can  readily  he  estimated,  .'■! ith  a knowledge  of  the  weighting 
function  (calculated  in  advance  by  ordinary  operational  methods) 
and  the  inpit  Xj.*  a suitable  time  increment  con  be  selected* 

The  arithmetic  operations  are  routine  snd  a standard  type  of 
graphical  computation  is  used* 

Having  obtained  on  approximate  solution  for  y(t),  it  Is 
still  necessary  to  find  ec(t).  In  some  caee3,  as  illustrated 
by  the  examples,  e0(t)  and  x2(t)  are  simply  related;  in  other 
cases,  a second  approximate  integration  may  be  required. 


Details  of  Solution 

Before-  presenting  any  exazaplos,  a sore  detailed  description 
will  be  given  of  tbs  method  used  to  obtain  the  response  of  the 
feedback  system* 

The  impulse  response  of  the  feedback  block  is  calculated 


by  taking 


/3(t)  = of  j>  (s)j  (20) 

Since  jS(e)  is  normally  the  ratio  of  two  polynomials  in  a,  the 
relative  orders  of  the  numerator  and  denominator  are  of  interest* 
If  yS(s)  is,  for  example,  of  the  form 


p > \ _ 

a \ » i — 


s " + 4 s + 15 


(21) 


the  Impulse  response  is  found  to  be 


^>(t)  = e sin  3t 


(22) 


from  any  transform  table0. 

If  the  numerator  and  denominator  are  of  the  same  order , as 


in  the  case 


/g(e)  = K -T 


(23) 


it  is  advisable  to  write  (by  division) 


/5u>  = 


3 + 1 


(24) 


The  impulse  response  is  then 


/5(i)  = K b (0)  - 
S (0)  is  a unit  impulse  at  zero  time* 


Eg 


(25) 


-_-p— -r  r-'Jj  k 


where 


I 

',3 
'«  - 

j 

I 


Such  Impulsive  components  actually  simplify  the  problem, 
since  they  imply  a 3trong  dependence  of  the  present  value  of 
~2(t)  on  the  present  value  of  v(t } . In  fact.  If  K S(0)  wore 
the  only  term  in  &(t)  t we  would  have 


Xo(t}  = 


y(t) 


(26) 


'Che  feedback  block  would  therefore  be  represented  by  a linear 
relation  between  the  Instantaneous  values  of  Xg  and  y,  and 
could  be  combined  with  the  nonlinear  forward  block  by  algebraic 
means. 

If  the  order  of  the  numerator  of  y2>(s)  exceeds  the  order 
of  the  denominator,  the  theorem  given  on  p<.  12  may  be  invoked » 
The  reciprocal  of  ^S(s)  is  then  placed  in  the  forward  block  and 
the  nonlinear  relation  is  inverted  and  put  in  the  feedback  block. 
This  possibility  is  Illustrated  in  one  of  the  examples. 

In  the  usual  case,  then,  a time  interval  h will  be  chosen 
and  the  values  of  it)  computed  for  t = k h,  where  k = 0,  1, 

2,  3*  . This  sequence  will  be  written 


A>.  P>>  fa . • • • 


(27) 


where  the  subscripts  correspond  to  the  various  values  of  k. 
Similarly,  the  values  of  y(t)  will  be  written,  as  the  sequence 


Jq9  ya>  yji 


(28) 


To  determine  the  values  of  x2(t)  at  the  times  t — n h, 
the  integral  of  2q.  (19)  will  be  evaluated  by  the  trapezoidal 


i 

i 

s 


$ 


rule,  Mcr?  elaborate  and  accurate  netho&s  are  arallable^"^ 
and  have  been  used,  in  similar  computations-'*  t but  the  choice 
of  the  trapezoidal  rule  results  in  a particularly  convenient 
and  simple  procedure*  According  to  this  rule , the  integral  is 


x2(llh)  ~ h I | A>  yn  * A yn-l  + A yn-2  + "* 


+ t^n-l  *1  * | /^n  -YoJ 

|^If  fy(t)  contains  an  Impulsive  component  of  magnitude  E,  Kyr 
should  be  added  to  Sql.  (29)  *1  Since  h jS^,  always  appears  In 
combination,  it  3 s convenient  to  absorb  h into  the  ^ -sequence 
given  in  (27).  The  first  member  of  the  /$ -sequence  can.  also 
absorb  the  factor  1/2 * as  can.  the  first  member  of  the  y -sequence, 
The  two  sequences  are  then 


Cf° . h/6l.  ^2.  h/b 


» yx  , 


(31) 


The  process  of  multiplications  and  additions  indicated  by 
Eq»  (29)  is  carried  cut  conveniently  on  a desk  calculator.  The 
sequence  (JO)  Is  computed  in  advance  and  arranged  in  tabular 
form  along  the  left-hand  edge  of  a piece  of  paper,  starting  at 
the  bottom.  The  values  of  y{nh)  are  determined  during  the  course 
of  the  computation  and  ere  recorded  in  a verliea]  -clurnn  on 
another  place  of  paper.  The  technique  used  to  keep  track  of  the 


mmUt 
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It  will  'oo  noted  that*  11111063  pQ  0,  .^{rfo.)  ‘will  depend 
ok  y , as  well  as  the  previous  values  of  y.  However , yr  -is  also 
dependent  on  %o > rh ; In  another  way,  as  indicated  by  trie  relations 

— •«*  rsV.l""!  {32i 


f ! x{nh)  s 

U J 


yn  f ^x(ch)  j (32) 

x(nh)  = x^Crih)  - x^lrsh)  (33) 

roia  Els-  (8)  e-sd  (9)*  If  we  denote  by  S(hh>  all  term©  of  Eq.  (29) 


except  -;ce  firot,  we  stay  write 


x,(r sh)  ss  ryw  4-  S(nh) 
^ n 


(34) 


where  r = h ^q/2.  Equation  {33i  may  therefore  be  written 


x(nh)  = x^(nh) 

- r y - S(nh) 

XX 

(35) 

= 3(nh) 

- r ?n 

(35) 

where 

3(nh)  = x1(nh) 

- S(nh)r 

(yn 

Since  x^Ct)  h&s  already  been  calculated  for  all  values  of  t and 
S Cilia}  depends  only  on  past  values  of  y,  th#  present  values  of  x 
and  y can  be  found  by  a sirnnlteneous  solution  of  Sis*  (32)  and 
(Ss)* 

The  simultaneous  determination  of  3t(hh)  and  y*,nh.;  can  be 
mode  very  conveniently  by  a graphical  method  shewn  in  Fig*  12* 
This  process  is  exactly  t-hs  ssme  as  the  one  used  for  the  large- 
signal  analysis  of  vacuum  tubes  having  resistive  loads*  The 
corresponding  quantities  are? 


2*  -g — — -?<-  :ty 

A* 


hammxtih  *3  the  slops  of  the  straight  line  given  by  2q«  (36)  is 
a c one t oat  and  only  d varies,  a single  straight  line  of  the 
proper  slope  can  be  dra-n  or  a piece  of  s®e! -transparent  tracing 
paper  which  may  be  superimposed  on  the  plot  of  y as  a function 
of  x.  As  ©oon  &e  §{ih!  has  been  computed,  the  tracing  paper  is 
placed  In  the  correct  position , the  values  of  y and  x(iii)  ere 
found  at  the  intersection  of  the  curves  and  entered  in  the  table 
of  values,  from  which  sf{n*l}hj  is  computed*  The  cycle  is  then 
repeated  as  many  tines  as  necessary « 


(I)  Series  Resonant  Circuit  with.  Nonlinear  Inductor 

A aeries  circuit  consisting  of  a resistor,  capacitor, 
and  nonlinear  inductor  is  shown  in  Fig.  13(a)  * The  problem 
is  to  determine  the  current  caused  by  a 5 -volt  step  input. 

The  original  block  diagraa,  given  in  Fig.  13(b),  is 
reduced  to  the  standard  form  given  in  Fig.  13(c)  by  moving 
the  block  containing  (l/e ) to  the  left  of  the  sutu-ing  point 
and  into  the  feedback  path.  Sine a 


S(s) 

— 

JL 

s 

(38) 

and 

^(a? 

— 

-S'  ■ 

(39) 

we  have 

X^(s) 

— 

(40) 

and  therefore 

^(t) 

= 

5t  . 

(41) 

The  feedback 

function 

is 

jS(s) 

= 

-2-  + -2z 
s 

(42) 

giving 

= 

2 -r  2t 

(43) 

For  b - 0.2,  the  modified  ^ -sequence  is  therefore 

, G.48,  0.56,  0.64  , 0.72  , 0.80  , 0.88,  **•  (44) 

The  quart if v r is  therefore  0»20,  and  the  general  equation  for 


?Ti  is 


;(ch)  - 0.20  y. 


"*ss*  mmmmmsnrn 
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To  check  the  accuracy  cf  the  method , a linear  csao  is 
created  first.  The  relation  between  td  and  i is  assumed  to  hs 

p i 

or  x “ y 


{46} 


I 

5 


The  initial  conditions  are  taken  to  he  q(o)  - 0,  1(0}  ~ 0* 

After  laying  out  the  table  sho  rn  on  the  next  page,  value# 
of  x,  '&?&  computed  from  Eq.  (41)  and  entered  in  the  x.  -column# 

* ri 

The  initial  value  of  jt  y0  = 0,  is  then  entered  at  the  top  of 
the  y-coluan.  The  next  entries  are  computed  as  follows? 

From  Sqs.  (29)  and  (44), 


S(0,2)  = 3X  yQ 

- (0.48) (0) 

= 0 

From  Eqs.  (37)  and  (4l) , 


S(0.2) 


^(0.2)  - S(0.2) 

1.000 


(4?) 


(48) 


Substituting  the  required  numbers  and  solving  Sqs.  (45)  and  (46) 
aitmlfcaneciusly , we  get 


a 


tg 


'4 

% 


I 


;! 


xCO.2)  - y (0.2)  ss  0.835 


(49) 


jj  The  graphical 

I this  uroblsa. 

) 

i 

| t-o  giv« 


cons  trust  ion  r shorn  in  Fig.  14, 
since  3qs.  (45)  (45)  can  be 


-Vi 


S(nh) 

1#2cT 


"*1 

is  unnecessary  in  I 
solved  analytic ally ! 


« 

i 


3( 


4 


i 


D 

0.400 

x.icrr 

2.036 

3.122 

4.304 

5,530 

6.792 

8.041 

9.259 


Table  I.  Calculations  for  Series 
with  Linear  inductor 


At  tlie  next  step,  we  li&va 
8(0.4)  s p\,  ; 


( 

I — 

3 ' 

*2 

x = ^ 

y = i 

1 eu 

0 

■ill.— 

0 

0 

1.000 

0.167 

0.833 

0.833 

1.600 

0.667 

1.333 

1.333 

1.893 

1.422 

1*576 

1.578 

1.962 

2.^65 

1.635 

1-635 

1.878 

3.435 

1.565 

1,565 

1.696 

4#59& 

1.414 

1.414 

1.452 

| 5.762 

1.218 

1 .213 

1.206 

6.993 

1.007 

1 »CZri 

0.959 

8.201 

0.799 

0.799 

0.741 

9.332 

0.616 

0.618 

Resonant  Circuit 


- Pi  J%  + Pz  /© 

- (0.46) (0.833)  + (0.56 )(G) 
» 

s'  0.400 


.-  ? -N  fc  % ...  . . \ **  i 6*  \ 

*/  * 2^ 


*"  2.000  - 0.400 

*■»  i 

x(0.4)  ~ y{0.4)  •-  1.333 


(51) 


(52) 

(53) 


jrr.v 

: r*  - 


M * 

if ' 


a *g  - 

r r * 
r *b  s 


2 *: 

n 

j 


f ^1*^1  * jr,-^!^ 


• *s  i 

tt  , 


At  the  next  step,  we  have 


s(o.6)  =t  8 x y*  + p2  y2 


p?>  y0 


= (0.48) (1.333)  « C0,56K0.i53»  + (C.64)(0) 


- 1.107 


(54) 


3(0.6)  = 


x» (0»6 ) 

A 

3. COO  ■ 
1*693 


biO.o; 


..lcr? 


x(0,6)  = y(C.6)  = 1.576 


(55) 

(55) 


*iiQ  rest  of  the  calculations  arc  made  in  similar  fashion, 
except  that  equations  like  (47),  (51),  and  (54)  are  handled  nj 
the  method  of  Fig.  11  rather  than  being  written  out.  explicitly. 

Similar  calculations  have  been  made  for-  h.  — O-i  and  h O.S, 
The  results  for  the  three  different  values  of  h are  given  in 
Table  II  on  the  following  page,  along  with  the  correct  values. 
Sfe©  j£«i  relation  assumed  is  such  that  the  inductor  has  an 
inductance  of  one  henry  for  the  assumed  single  turn.  The 
r-elatlon  between  the  current  aud  applied  voltage  is  easily 


shown  to  be 


1 1 3) 


For  a 5-volt  step,  it  follows  that 


2s  + 2 


(57) 


i(t)  = 5 © sin  t 


(56; 


/*£“**" 
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Current 

f 

Approximate  Values 

| Correct  Value 



... 

fc.  “ 0.1 

_ _} 

h “ 0 ,2 

j n - 0.5 1 

0 

0 

0 

0 

F 

0 

0.1 

0.454 

s 

: 

0.452 

0.2 

0.816 

j 

0.834 

0.614 

o.i 

1.101 

1.095 

0.4 

1013 

1.353 

1*>T: 

0.5 

1.451 

1.567 

1 ,*0-4 

0,6 

1.557 

1.578 

3 .550 

0.7 

1:606 

1.-500 

0.8 

X.61S 

1.635 

1.612 

0.9 

1.559 

1.593 

1,0 

1.552 

1.565 

1.567 

1.543 

, o 

A •- 

1.414 

1,404 

1.4 

1,218 

1.215 

1.5 

1.111 

1.113 

X .6 

■*  r,iV.' 

1 *W| 

1.009 

1.8 

0.799 

0.805 

2.0 

0.616 

I 

i 

0.556 

0,615 

. 

o C 
- o 

1 

0*186 

0.246 

}.o 

‘ 

i 

0 000 

0 036 

■ 

v c 

^ . -> 

-0.061 

-0.053 

- 

4.0 

( 

—O.G62 

-C.C69 

1 

-0.042 

-0.055 

! 

■ 

1 

5.0 

| 

-0.022  | -0.032 

Tzhlst  IX . Hebuita  of  Calculations  for  Series  Besoz 
Circuit  JIth  Linear  Inductor 
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■4  i>Iot  of  the  values  in  T&oxe  11  shows  that  the  curve 
for  h “ v*".  is  practically  coincident-  with  the  correct  curve. 
Tb-i  - " ;«fe  'i-.T  h ~ 0.2  differ  ’ey  not  wore  tnan  0.030  from  the 


.S\  vrte  ran*;s 


an  error  of  not  more  than 


.•  .-a'.  : per-  C v XT'S*  4 i£2£5?22.' v OX  .CL  0*3  it»ads  to  fairly  big 

i-r  but  gives  a reasonably  good  qualitative  pictures  the 
- '.ctirasita  and  minimum  points  of  the  curves  occur-  at  the  proper 

-V 

; .usees . In  general  t the  error  varies  roughly  as  hX  * 

The  same  process  is  employed  if  the  inductor  is  nonlinear, 
except  that  y and  x are  no  longer  related  by  a simple  pro- 
portion* For  the  example  illustrated  here,  the  curve  of  Fig.  15 
is  used*  Remembering  that  y = i and  x = ps  it  will  be  seen  that 
this  curve  is  a typical  magnetization  curve , applicable  to  a 
variety  o?  common  magnetic  materials. 

Suppose  that  h ~ 0*2  ia  taken  for  the  time  increment.  The 
same  table  is  used  for  the  ,6-sequence  and  a table  lihe  Table 
11  is  started*  beginning  with 

3(0.2)  = X.OOQ  (39) 

as  before.  In  this  case,  as  indicated  by  Fig,  13,  the  values 
of  x and  y ar  c 

x(0*2)  ss  0-89 


y(0»2)  = 0.57 


(6o) 


If  values  of  current  (y)  are  gill 
no  need  to  record  x,  since  it  is 


that  is  required , there  is 
not  used  in  the  computations.^ 


WIMUHI4rMtVlllh  iO**-” 


^ id 

•'*  S-J  ** 

At  ti;o  ne: 

si  step,  the  values  ere 

3(0.4) 

--  (0.48)  (0*57)  + (0.55  )(0) 

- 0c2?4 

(61/ 

3(0.4) 

= 2.000  - 0,274 

h 

= 1.726 

(62) 

x(0.4) 

y(0.45 

--  1 . — 0 

(from  Fig.  15) 

- 2.35 

(63) 

d i 


The  calculations  are  continued  In  similar  fashion  as  far  as 
necessary*  Results  of  such  calculations  for  h = 0.05,  0.1  and 
0.2  are  given  in  Table  III  on  the  following  page  and  plottel 
In  Fig.  16. 

Inspection  of  the  table  and  the  curves  indicates  that 
h = 0.2  is  not  an  unreasonable  increment,  except  for  one  poor 
value  at  t = 0.4.  For  small  values  of  t,  the  small  Increments 
give  better  results.  In  all  cases  there  is  some  fluctuation  in 
the  values,  caused  In  port  by  inability  to  obtain  more  than 
throe  significant,  figures  from  Fig.  15.  Some  of  the  accuracy 
attainable  by  the  use  of  the  small  XTiCreiaei.ua  i3  not  realized 
for  this  reason. 
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O 

OJ 

n,2 

0.3 

0*4 

0.45 

0.5 

0.6 

0.7 

0.8 

0.9 

1.0 

1.1 

1.2 

1.3 


• 9 

X s*t 


XO 

1.6 

1.7 

1.8 
1.9 


o n 

— ■ • ■©  v 


i 

i 

Approx 

I .n  ate 

Current  V e', 

1_U  3 = 

L 

? 

i 

h = C 

>.05 

H 

i 

11 

o 

• 

M 

i 

L 

. - 0.2 



t 

i 

I 

f'K 

A 

v 

0 

\ 

j 

i 

n.10 

0,10 

f 

\ 

0.69 

i 

0.54 

0.57 

1.67 

1 

; 

• 

1.80 

i 

2.06 

2.10 

2.05 

1 

! 

i 

2 .11 
2*06 

2.35 

1.90 

1.88 

1.90 

1.74 

1.74 

1.55 

1.56 

1.55 

1.38 

1.25 

1,22 

1.12 

0.93 

0.97 

0,90 

1 

| 

0.80 

0.80 

» 

5 

| 

1 

| 

0.73 

1 . 
i 

i 

\ 

i 

s 

0.65 

0.60 

0.65 

j 

i 

0.54 

0.54 

V 

i 

3 

1 

I 

i 

! 

* 

1 

\ 

r 

? 

0.48 

0.44 

1 

( 

l 

1 

0,43 

fable  XII.  Results  cf  Calculations  for  Series 
Resonant  Circuit  with  Konlineax*  Inductor 
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(2)  PC*  kafijji'-  Heivrorb  Containing  a Sleds 

The  circuit  diagraia  of  the  SC  ladder  network  considered 
in  the  sscrjiiv  example  is  sho.-m  xu  Pjqg,  17(a);  the  block  diagram^ 
is  given  in  Pig.  17(b) * 


1 

\ A A . ..  .... 

V V V j 


_L 

ci  ~p 


1 


^ «2 

L, 


Figure  1T»  HC  Ladder  network  sssd  3.lock  Diagr.ym 


In  this  oas«p  the  diode  Is  described  by  a conductance  function, 

•ii 

denoted  by  ^ • The  variable  y la  the  current  'throu^ti  the  diode* 


1 

e i 

r~ — 

X 

k 

R-jCtS  + 1 
1 * 

i 

i R 


I B3 

'pi  R2C,2s  + 1 


I I 

% i ! 


PlCIs  + 1 %C2®  +1j 
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For  simplicity  , it  is  assumed  that  R,  ™ Rp  = 1 and  that 
0^  = 0^  = A 3-voli  step  is  applied  at  the  inpjt  terminals 
«a  the  network  with  the  capacitors  initially  uncharged* 

'-nth  these  values , the  preliminary  calculations  are  as 
follows*  The  transform  of  the  Input,  voltage  la 


e(s)  = — 2™ 

e 

and  the  transfer  function  of  the  inpat  block  Is 


with  tne  result  that 


G^s) 


x-.(s)  = 


8 + 1 


s (s  + 1) 


Xx(t)  = 3 (1  - e } 


The  transfer  function  cf  the  feedback  block  is 

2 


(s)  - 


s + 1 


so  that 


d 3 


^(t)  - 

The  is  therefore 

^32,  0.2562,  0.2196.  0.1798,  *'*• 


(64) 


(65) 


(66) 

(67) 


(63) 


(69) 


(70) 


if  h ss  0.2(5® 

It  may  also  be  noted  that,  for  the  particular  values 
selected. 


&aa  therefore 


<x  _ xa(s) 
o 2 

(71) 

Xb(t) 

e (t|  = — *s—  * 

(72) 

I 


i trw»rXMiSr*=^ 


’ f4*AiA^  'f  •,-■,-■*■  >«.,  «• , 


A linear  case  is  censi dared  first , as  a cheek  on  the 
accuracy  of  the  method * If  a 1-ohm  resistor  is  substituted 
for  the  diode 9 

x = y (73) 

Proceeding  as  before,  fable  IV  is  obtained. 


i 

i 

t j 

X*s 

3 

s 

■ 

x2 

e 

0 

A 

1 

r — rrr“^ 
0 

1 

- 

- 

0 

0 | 

0 

0.2 

j .5438 

.0000 

.5438 

.4532 

.0906  i 

1 

.0453 

0.4 

| .9890 

.1484 

. 8406 

.7005 

.2885  ! 

j 

.1442 

0.6 

j 

i 1.3536 

i 

<>3509 

1.0027 

.8356 

.5180  1 

.2590 

0.8  1 
i 

1.6520 

.5610 

1.0910 

•9092 

.7428  : 

.3714 

1.0 

1.8964 

.7571 

1.1393 

.9494 

.9470 

i 

.4735 

1.2 

2.0964 

.9308 

1.1656 

• 9713 

1.1251  ] 

.5625 

X»4 

2.2602 

1.0602 

1.1800 

.9833 

1.2769  I 

.6  >54 

1.6 

2.3943  | 

1.2064 

1.1879 

.9899 

1.4044 

.7022 

1.8 

j 

2.5041 j 

1 .3118 

1.1923 

.9936 

1.5105 

1 .7553 

2.0 

1 2.5940  | 

1.3993 

1.1947 

.9356 

1.5984  | 

i .7992 

1 

Table  IV.  Calculations  for  RC  ladder  f et'ork 
with  Linear  Resistor 


Similar  calculations  were  made  for  h — 0.1  end  h = 0.5* 
The  results  of  three  series  of  calculations  are  summarized 
in  Table  V,  'which,  also  includes  the  correct  values  computed 
from  the*  relation 


9 { t)  = 1 


e * + 


e 


>t 


(74) 


».«gg"  «*■ 
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n. 


Output  Voltage 


t 

* 

tb±  i U* 

■•exinate  Values 

Correct  Value 

i 

i 

1 

h = 0.1  [ h = 0.2 

| h - 0.5  “ 

0 

0 i 

1 

! 

0 

r 

r = ~ 

! 0 

0.1 

0 . m ■Eft  f 

! ~ — l 

0.0132 

o.a 

0.0461  ! 

i ! 

0.0453 

0.0463 

0.3 

| 0.0918  i 

1 

0.0921 

0.4 

0.144-9  | 0.1444 

0.1451 

G»5 

0.2036 

c.1967 

0.2018 

0.6 

0.2594 

0.2590 

I 

i 

1 

0.2594 

0.7 

0.3163 

0.3164 

0.8 

0.3713 

0.3714 

(\3714 

0.9 

0.4238 

0 - 4 2 3 2 

1.0 

0.4732 

0.4735 

0.475c 

0.4731 

1.2 

0.5525 

0.5619 

1.4 

0.6384 

0.6376 

1.5 

0.6762 

0.6709 

1.6 

0.7022 

0.7013 

1.8 

0.7553 

0.7543 

2.0 

0.7992 

0.8046 

0.7982 

3.5 

0.8847 

0.8771 

3.0 

! 

i 

1 

! 

0.9329 

0.9253 

3.5 

! 

0.9617 

0.954? 

a# 


Table  V'.  Results  of  Calculations  for  30  Ladder 
Netv/ork  v;lth  Linear  Resistor 
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An  an  approximate  on  to  the  volt^;e~current  characteristic 
of  a crystal  diode,  the  relation 

1 = 0 

4 

& 

was  used.  In  the  general  notation  used  previously, 

* 

y = x ' x > 0 (76 } 

Equation  (76)  was  plotted  and  the  calculations  were  repeated, 
using  the  graphical-  method  «*_.  .rai’act  before  Results  for  h — 
0.1,  0.2,  and  0.5  are  given  in  Table  VI  and  plotted  in  Fig.  16. 
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t 

Approximate  Output  Voltage 

h = 0.1 

h = 0.2 

h = 0.5 

0 

0 

0 

0 

0.1 

0.000, 

i 

7 

| 

0.2 

0.005 

0.008 

0.3 

0*022 

0.4 

0,052 

0.068 

0.5 

0.121 

0.150 

0.6 

0.190 

0.192 

C.7 

0.261 

! 

0.8 

0,329 

C.?28 

0.9 

0.392 

1.0 

0,450 

0.450 

0.442 

2.0 

0.800 

1 

0.800 

e < 0 
e > 0 


(75] 


Table  VI*.  Results  of  Calculations  for  RC 
Ladder  Network  with  Nonlinear  Diode 


■«75^sr 


bO 
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In  the  linear  case,  the  solution  for  h =s  0.1  differ: 
from  the  correct  solution  by  not  more  than,  three  units-,  i 
the  fourth  decimal  piece.  For  h = 0.2  and  Q.5>  the  msxXss 
errors  are  10  and  80  units,  respectively,  As  before,  the 
error  “arias  roughly  as  h2  - 

Since  a correct  solution  c-o  the  nonlinear  prx-Lieai  i 
not  available,  the  error  c rennet  ha  determined  preoisely* 
however,  it  apoeara  that  tbs-  results  for  h = 0.1  are  sbe- 
as  good  sb  can  b-:  expected  with  the  number  of  significass: 
figures  enspleyed. 

(3)  Torque -Saturated  Servoaschanl.esi 

As  another  example,  the  simple  torque -ssuur &tsc  « - <* 

mechanism  shown  in  Fig.  19  is  consider 


s(Js  t f ; 


Figure  19 « Block  Oiagram  of  Torque-Saturated  uervomechtinieis 


ha 


The  block  diagram  of  the  torque -saturated  servomechanism 


is  reduced  to  the  general  form  oj  moving  the  linear  block  to 
the  right  of  the  take-off  point  ana  into  the  feedback  path; 
the  transfer  functions  are  then 


Ma) 

A 

(s)  = G0(s> 


5 (J3  + f) 

The  relation  between  error  and  torque  may  be  expressed  as 


in) 

(78) 


= ktc 

= X £ 


- K £. 


e>£c 

ec<e  <ec 

€<-ec 


(79> 


Without  loss  of  generality,  a new  set  of  dimensionless 

22 


parameters  may  be  Introduced;  a suitable  choice  is 

f 


zrfT: k 


(80) 


05  _ — 


~€  = e«0  t 


/ ao  * 
\v*~  / 


(83) 


In  the  example  which  follows,  w©  have  chossn  to  let  S ~ 0.5, 
= 1,  and  Ce  - l.  This  is  equivalent  to  selecting  K — 3 f 


&> 


;T  = 1,  «r>d  f = 1 in  the  original  system,  and  leads  to 


% = 


x 


e 


(84) 


(85) 


iLhf.  „ 

It  will  be  assumed  that  a 5-unit  step  is  applied  tc  the 
input  with  the  system  initially  at  rest.  For  the  values  used, 
w©  havo 


xi(t> 


/8 


(*  'i 


- 1 
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X j 
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Si?; 
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for  h — 0.4,  the  ^3 -sequence  becomes 


(88) 


0.0000,  0.1319,  0.2203.  0.2795,  0,3192,  *••• 


(89) 


The  necessity  for  the  graphical  solution  of  two  algebraic 
equations  is  eliminated,  since  & (0)  “0.  In  this  particular 
situation,  the  value  of  y at-  any  time  is  dependent  only  on  the 
past  values  of  x and  can  be  calculated  without  a knowledge  of 
the  present  value  of  y.  The  calculations  are  started  as  follows; 


At  t«  — Qt 
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x(Q) 


5.0000 

0.0000 

5=0000 
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*2(.4}  = (6.1319)  Uj£G&2) 

ss  0.0660 


%{ *4) 
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At  t s 0.4; 
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Theae  calculation®  are  best  carried  out  in  tabular  form? 
part  of  the  work  is  given  in  Table  VII.  As  a reminder  that 
y(0)/2  ia  used  in  the  calculation  of  x2(t)»  the  correct  value 
of  y at  % & 0 has  been  crossed  cut-  am  one -half  of  the  correct 


value  entered  at  this  point  in  the  taoie. 
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k2 
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y = T 

0 

0 
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5.0000 
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b.oooo 

XirWW 

0.4 

(1 

0.0660 
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4-S340 

1,0000 

0.8 

1.2 

« 

15 

0.2421 
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1*6 
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0.7908 

4.2092 
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2,0 
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it 
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1.1228 

3,8772 

1.0000 

1 

2,4 

1,4776 

3.5224 

l.OOCO 

2.8 

?? 

1,8471  | 

3.1529 
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2,2265 

2.7731 

l.*OOOC 
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u 

2.6132 

2.  >866 

1.0000 

4,0 
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3.0042 

1.9958 
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4.4 
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3 . 3961  j 

1.6019 
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1. 0000 
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1..2G3C 

1.0000 

5.2 

f* 

j 4,1912 

1 
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Table  VII,  Calculations  for  Torque-Maturated 
3 ervcmaeh&ni  am 


calculations  have  bean 


f j SKiSM 
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ff* 
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lit  Fige  £0,  the  results  of  the  calculations  have  bean 
superimposed  on.  a <nar*re  which  was  determined  analytically*^  * 
For*  t <6,  the  approximate  values  obtained  by  the  step-by-atap 
method  are  consistently  higher  than  the  c err act  values | the 
difference  is  not-  more  than  140  units  in  the  fourth  decimal 
piece*  In  a.  separate  set  of  calculations f not  included  here, 
h js  0.2  was  used?  the  difference  between  the  approximate  and 
correct  values  does  not  exceed  40  units.  For  t >6,  the 
accuracy  is  somewhat  better. 

A series  of  calculations  was  also  made  for  en  error- 
torque  relation  which  may  bo  expressed  as 


and 


r 

K%  i1 


t(-£)  = - -x(e  ) 
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-(4-)' 

*-'C 


S>  0 


!g$) 


i i,*4  j 
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This  expression  for  the  torque  reduces  to  X £ for  small 
error’s  am  approaches  K £0  for  large  errors,  life©  the 
idealised  51  linear"  or  "straight-line"  saturation  characteristic 
used  in  the  earlier  calculations . Th.a  resulting  curve  for  the 
cavui  avion  cn.noc tcristxc  ic  wuussc  plotted"  in 
Fig.  2o.  mere  is  very  little  difference  in  ana  curves  for 
iax^e  ezTora > in  the  neighborhood  of  zero  error,  the  linear 
saturation  characteristic  leads  to  a larger  overshoot,  as 
might  be  expected . 
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V“t  » rsO  KstworiC  t:  oi'iw  arguing  s.  iJioo-C 

The  final  sxa&pie  in  this  report  i®  the  EC  network  shown 
in  Fig.  21<&).  If  the  block  diegraa  is  contracted  and  reduced 
to  the  general  fora.  Fig*  21(h)  is  obtained*  Since  /5(&)  has 
a numerator  of  higher  order  then  th-s  denominator  , it-  will  be 
advantageous  to  interchange  tho  forward  and  feedback  blocks 
by  use  of  the  theorem  descri  bed  os  psg©  jk  j the  resulting 
block  diagram  Is  given  in  Fig.  21(c)* 

She  step-by«3tsp  calculation  process  for  systems  having 
a block  diagram  of  the  type  given  in  Fig*  21(c)  is  similar  to 
the  process  already  described,  but  cLi?  f in  some  details. 

A general  discussion  of  the  required  modifications  in  the 
method  will  therefore  be  presented  before  this  specific  case 
is-  treated* 

For  a feedback  system  with  the  nonlinear  block  in  the 

feedback  path,  the  basic  equations  are 


3t-,(nn)  ■=  x,  (hhj  - x(nh) 

CL  1 


x(rh)  = 


r , .1 

g £ y(nh) i 

i-t  ««K 


y<nh)  = j x?(t  -t)  a(X)  atr 
s ~r*^  x-.  (nh)  t Q(nh 


(95) 


(96) 
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•98} 
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(99) 


and  Q(hh) » like  the  S(nh)  usod  previously , denotes  a weighted 
nMn  of  past  values  of  , up  to  [(a-X)iil  * obtained  as  part 
of  the  approximation  to  the  convolution  integral. 
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ombinlng  Eqs,  (95)  and  (98)  ard  sotting  Ip  /2 


*»t-  obtain 
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yCsh)  = p I x<|  (r4*)  - x<na)  i 4-  qChh) 
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where  4 (nil)  ~ H(nh) 


A>  x-  (rii' 


002} 


Equations  (96)  and  (101)  constitute  a pair  cl  simultaneous 
algebraic  equations  in  x and  y which  osn  be  solved  by  a 
graphical  process  illustrated  in  Fig;,  22(«) » Alternatively, 
to  permit  U3e  of  a plot  of  y as  & function  of  s,  the  axes 
car?  bt?  interchanged  as  indicated  in  Fig-  22(b). 

For  this  example ; it  will  be  assumed  that  Ri  ~ 1,  C1  = X, 
and  the  input  is  a 2 -volt  step,  applied  with  the  capacitor 
initially  uncharged,  A linear  case  for  which  x = y,  obtained 
by  replacing  tho  diode  by  a linear  2 -cbm  rooistor,  will  be 
treated  first.  With,  these  particular  values  and  using  h = 0,2, 
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At  b = 0,20?  Q( 0,2)  = 


r.  0\ 

V i- . 


- 0.163? 

= Q(0.2)  + <t  x*  ( 0 .2 } 

f 

t 

=s  0.1637  * 0*2000 

0.363? 

Using  the  graph * x(0.2)  - y(0*2)  - 0.3306 
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us: 


_ _ _ i 2.  u « / r V i,  V 

,iis?  ncratsi-  s;  trt  y.4  j 


~ y «.  u ; 

x^(0.4)  ~ 2*0000 


0,5522 


0,5522 


1,4476 

Lot's  are  carried  out  in  similar 
fashion;  the  «o£*k  is  summarised  in  Sable  VXXX*  In  this  tsKsa, 
&s  in  fable  VII , the  value  of  y{€\5  is  crossed  out  and  replaced 
by  y(0)/2  to  facilitate  the  subsequent  mult5.pl  is  at  Ions  which 
determine  the  values  e£  Q,* 


VIII . Calculations  for  BO  K«tw*or!r  with  Linesi 
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The  errors  e.t  t.  ~ 0*8  are  6 
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fourth  decimal  place?  the  variation  is  again  roughly  as  hf  * 
in  this  case,  the  values  obtained  using  *-.  0*4  appear  to  he 

quits  satisfactory  for  plotting  purposes, 

A similar  sot  of  calculations  was  made  fir"  this  network 
with  a nonlinear  diode*  described  by  the  relation 
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1.3.0  J 


Sines  y is  identified  with  output  voltage 
«q*  (113)  becomes 

jm  j—m  are-"'* 


id  x with  current. 


(114) 


in  the  stersArwi  notation.  Following  the  scheme  of  Fig*  22(b), 
the  graphical  solution  looks  exactly  the  same  as  in  Fig.  15, 
except  that  x and  y are  Interchanged*  The  results  are  gives 
in  Table  X. 

In  this  problem,  the  calculations  with  h = 0®4  appear 
to  exhaust  the  accuracy  of  the  graphical  solution.  Use  of 
smaller  steps  doss  not  increase  the  accuracy,  since  the  main 
source  of  error  is  the  graphical  determination  of  x and  y* 
The  smaller  steps  do,  of  course,  give  & more  complete 
picture  of  the  response. 

The  results  are  also  plotted  in  Fig.  23,  using  stops 
of  h =s  0.1  for  t<l*0  and  h ~ 0*2  for-  t>1.0*  The  correct 
(computed)  curve  for  the  linear  case  is  included  for  com- 
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The  accuracy  ox  the  results  obtained  by  the  procedures 
outlined  in  this  reoo.rt  is  affected  bv  c.  variety  cf  possible 


errors : 


(1)  O&iculatica  errors  - due  to  arithmetic  slat ekes 
in  carry  ing  out  the  various  r-ult  iplic at ions , 
additions,  and  subtractions 

(2)  <>rar4xlcal  errors  - due  to  inability  to  specify 
the  intersection  of  two  curves  more  closely  than 
three  or  four  significant  figures 

(3)  Round-off  error 

(4)  Integration  error  - due  to  inadequacies  of  the 
approximate  integration  formula  (trapezoidal 
rule) t assuming  that  the  values  of  the  integrand 
are  correct 

(5)  Inherited  error  - due  to  inaccuracies  in  the 
integrand  resulting  from  previous  errors 


The  possibility  of  errors  of  the  first  three  types  can 
be  minimized  by  elementary  stratagems  s (1)  careful  work  and 
different; e checks  of  the  variables  as  they  are  obtained;  (2) 
use  of  a sharp  pencil  and  a large  pise©  of  graph  paper;  and 
(3)  carrying  mere  decimal  places  than  will  be  retained  in 
the  final  results.  The  other  types  of:  error  deserve  a more 
detailed  discussion* 

Bsc s&ise  the  trapezoidal  rule  is  based  on  an  approximation 
of  the  integrand  by  a series  of  straight  lire  segments,  it 
will  give-  correct  results  for  any  integrand  having  a constant 
first  derivative;  any  errors  which  are  encountered  can  thus 
fee  traced  to  the  higher  derivatives  ox  the  integrand.  A 
thorough  analysis  of  the  question  leads  to  tho  following 
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estimate  of  the  error? 
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(115) 


v/fceye  fe  sc  size  of  increment  between  values  j 

L ~ iih  ~ tots?  interval,  of  integration; 

s s function  to  be  integrated  (integrand); 

f'JL  = maximum  value  of  the  second  derivative 

of  the  integrand  in  tha  interval  0<t<nh. 

A more  optimistic  estimate  of  the  error' 


[- 


(rih)  - z'CO)  (116) 


where  z. ? ( m ) and  z'  (0)  are  the  first  derivatives  of  the 
integrand  at  the  end  and  beginning  of  the  interval*  Since 


the  first  derivative  is  the  integral  of  the  second 8 it  is 
true  that 

rnh 

z”  dt 

J Q 

= 1 Z"m  (117) 

wh  are  z „„  denotes  the  sverege  value  cf  the  second  derivative 

E? 

of  the  integrand  in  the  interval  under  consideration#  The 
two  formulas  ere  therefore  essentially  the  same.  The  m±mz?> 
sign  in  the  formulas  reflects  the  fact  that  the  eipproxismte 
value  of  the  .integral  will  bo  too  high  if  the  curve  which 
represents  the  integrand  is  concave  upward,  placing  the 
straight-line  approximation  above  the  curve# 
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‘fhe  integral  to  be  evaluated  la  either 
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or 
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depending  on  the  nature  of  the  block  diagram.  If,  for  ©x«Eples 
we  define 

zCt)  = &(t)  y(t  -t)  , (118) 


wo  then,  have 


z*it)  = 


*{t)  y*( t 


«v  % 


+ «*  it)  y(t  - 'l'  ) 


and 


«’s  ' 


(119) 

fi  00  '/" (t  -T)  - 2 <6*00  7 s (t  -U) 


+ ^8“  00  y(t  -TC)  , 


(120) 


where  the  prices  denote  differentiation  with  respect  to  T . 
While  j8(t)  or  •/ ( t ) is  known  in  advance*  y(t)  and  Xg(t)  are 
determined  in  the  course  of  the  calculations  and  are  subject 
to  the  orrora  hied  are  to  be  estimated.  day  estimato  of  the 
error  must  therefore  be  made  after  the  calculations  have  been 
completed.  Equation  (119) » in  particular » car  bo  usod  to  obtain 
a quick  estimate  of  the  error  as  the  calculations  progress. 

Perhaps  the  must  useful  ecncluaic n to  be  drawn  from  fiqa. 
(115)  and  (1 26  i Is  the  variation  of  error  vr ith  h.  If  two  series 


of  values  have  been  for  s -3! ven  v-ar?  able  by  using  tv/© 


different  valued  of  hs  and  if  these  values  are  net  subject  to 
other  errors  f the  fart  that  the  error  varies  as  hr  night  bo 
used  to  pick  a third,  value  of  h which  would  lead  tc  results 
of  any  desired  accuracy. 
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Although  the  question  of  inherit el  errors  has  not  been 
thoroughly  studied,  the  following  observations  ear:  be  made. 
Suppose  that  the  system  is  described  by  the  relations 
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y(t>  = f 


.ail 
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The  value  of  Xg(t)  io  approximated  by  a v ’sighted  sum  of  present 
and  past  values  cf  y.  If  the  approximate  integration  process 
introduced  no  error  except  that  due  to  inaccurate  values  of  y% 
the  error-  in  would  be  a weighted  sum  of  the  errors  in  y. 
Assuming  that  f (>:)  la  a well-behaved  function  and  that  h is 
small  enough  so  that  there  aro  only  small  changes  in  :Cj(t ) and 
|S(t)  in  any  increment , there  is  a tendency  for  errors  in  x 
and  y to  alternate  in  sign.  This  tendency  is  evident  in  some 
of  the  examples  and  may  be  explained  as  follows.  Suppose  that 
one  of  the  values  of  y Is  too  large j the  corresponding  value 
of  Xg  from  3q„  (19)  will  likewise  be  too  large.  The  next  value 
cf  x,  computed  from  3q»  (8)^  rill  bo  too  small  ard  will  lead 
to  a value  cf  yf  from  Tc , (9),  which  will  therefore  bo  small , 

If  the  errors  are  not  too  big,  the  process  thus  tends  to  be 
partially  self -correcting,  with  the  result  that  inherited 
errors  have  not  been  a problem. 

If  the  step-by-step  process  is  to  be  applied  to  physical 
systems  and  checked  experimentally  .,  other  sources  of  error  aro 
Introduced;  the  most  important  is  th-a  possibility  that  dynamic 
behavior  of  nonlinear  elementa  is  not  predictable  from  the 
static  characteristics  • 


w 
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In  addition  to  knowing  about  the  aecu rsoy  of  & step- 
by-step  mothod,  the  user  also  want  a an  estimate  ex'  the  amount 
of  computational  labor  requi red  to  c caplet®  & solution,  In  the 
method  described  here,  the  Ic®1  step  in  the  calculation  will 
require  about  (2  + k)  operations,  either  srultiplic at  lens , 
additions,  or  subtractions,  and  one  dot ercinst  ion  of  values 
from  & curve  or  equation*  Summing  o^er  n steps,  the  number 

of  operations  is  e» l- lasted  to  be  approximately 

n (n  + 1} 

3 n + ■ '-»•■ » -»  - —■ » ■ - ~ 


cm) 


which  does  not  include  the  preliaJjiary  calculations  needed 
to  obtain  x (i)  and  &Ct).  Since  t = hh,  the  estimate  may 

also  fe®  expressed  as 

- 2 «• 

Ail  'y  '■  A i™  A'——  ( 11^ } 

2 h"  2 fc 

indicating  that  the  errors  (which  vary  as  h2 } will  be  an 
inverse  function  of  the  effort  expended. 

This  ©c-tlmst©  compares  r&t her  unfavorably  with  ten.®*5 
approximate  methods  for  which  the  operations  per  step  is  a 
constant.  In  defense  of  the  proposed  method,  it  may  fee  said 
(I)  that  the  ramiber  of  operations  is  inicpeii&ent  of  the 
order  of  the  aquation  being  solved  and  (2)  that  the  use  of 
the  entire  past  history  of  one  of  the  varia&laa  helps  to 
ssfca  possible  the  sccur-sciea  attainable  by  this  method.  In 
special  caaoa,  if  /3(t)  approaches  a constant  for  large  t, 
soma  of  this  work  may  tea  avoided. 


-cm'  SS^X*V:*t^*e?:£,S 
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A method  for  the  approximate  alternation  of  the  response 
ef  a nonlinear  system  to  fairly  arbitrary  inputs  has  been 
described.  The  method  is  presently  applicable  to  systems  with 
©so  essential  nonlinearity  whose  mathematical  relations  car-,  be 


indicated  by  & standard  type  of  block  diagram-  The  techniques 
involved  — operational  mathematics,  rumor  leal  integration,  and 
graphical  solution  of  algebraic  equations  — ars  generally 
known  to  engineers,  who  will  find  the  accuracy  of  the  method 
satisfactory  for  most  purposes* 

It  appears  that  the  method  is  applicable  to  systems  in 
which  the  nonlinearity  changes  in  a prescribed  manner  with 
time  or  depends  on  an  independent  signal  * such  as  a grid 
voltage  of  a vacuum  tubs,  tflth  some  additional  work,  it  may 
prove  applicable  to  - systems  having  more  than  one  nonlinear 
a lenient. 
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